Mon. Not. R. Astron. Soc. 000, [HS] (2010) Printed 13 July 2011 (MN WT^ style file v2.2) 



First observations of CN(2-1), HCO+(3-2) and C2H(3-2) 
emission lines in the Perseus cluster: constraints on 
heating mechanisms in the cluster gas. 
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ABSTRACT 

We present the first observations of emission lines of CN(2-1), HCO+(3-2) and C2H(3- 
2) in the Perseus cluster. We observed at two positions: directly at the central galaxy, 
NGC 1275 and also at a position about 20" to the east where associated filamentary 
structure has been shown to have strong CO emission. Clear detections in CN and 
HCO"*" transitions and a weak detection of the C2H transition were made towards 
NGC 1275, while weak detections of CN and HCO''" were made towards the eastern 
filamentary structure. Crude estimates of the column densities and fractional abun- 
dances (mostly upper limits) as functions of an unknown rotational temperature were 
made to both sources. These observational data were compared with the outputs of 
thermal/chemical models previously published by iBavet et al.l ()2011[ ) in an attempt 
to constrain the heating mechanisms in cluster gas. We find that models in which 
heating is dominated by cosmic rays can ac count fo r the molecular observations. This 
conclusion is consistent with that of Fcrlan d et al.l (|200 9) in their study of gas traced 
by optical and infrared radiation. The cosmic ray heating rate in the regions probed 
by molecular emissions is required to be at least two orders of magnitude larger than 
that in the Milky Way. 

Key words: Astrochemistry - ISM:abundances - galaxies: intergalactic medium - 
Galaxies: clusters: individual: Perseus - galaxies: individual: NGC 1275 - Galaxies: 
cooling flows 



1 INTRODUCTION 

The Perseus Cluster is one of the nearest galaxy clus- 
ters and is the brightest X-ray cluster in the sky. 
The cluster and its central galaxy NGC 1275 have 
been the focus of intense study for many years, at 
X-ray, optical, IR and millimetre wavelengths. The 
first molecular detections were of CO ro tational emis- 
sion toward s the centre of NGC 1 275 JL azarcff et al.! 
1989; Mirabel. Sanders fc Ka^ Il989l: fReuter ct al. 1993,; 



ture observed in Ha |Hu et al. _ 19831 ; ICowie et al.l 119831 ; 
IConselice. Gallagher fc Wvsell200l|) within the hot gas de- 



Braine et al.lil995l;llnoue et al.ll 1996 1; [Bridges fc Irwinlll99i ; 
Lim. Ao fc Dinh-V-Trung I2OO8I I. Observations have now 



demonstrated that CO emissio n also extends to some te ns 
of kp c from the central galaxy (jSalome et al.ll2006l . l2008al lbl. 
boilh and is strongly correlated with the filamentary struc- 
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tected in X-rays at 0.5 keV (|Fabian et al.l HoOS). The 
Ha structure is also correlated with warm H2 emission 
lEdge et al.ll2002l ; IWilman et "alll2002l '). Some of the molec- 
ular ga s towards the c luster must be at high density (^ 10^ 
cm~^'l. ISalome et al.l (|2008al ) have made the first detection 
of the high density tracer HCN(3-2) emission towards the 
centre of NGC 1275. 

The nature of the filamentary s t ructur e is the subject 
of current discussion. I Salome et al.l ([20111) note two pos- 
sibilities: either that the CO filaments form far from the 
galax y's centre from upl ifted warm gas, eventually falling 
back (jRevaz. Combes fc S alome 2008), or that the molecu- 
lar gas is entrained and dragged out of the galaxy by ris- 
ing hot gas. Evidently, further observations of molecular gas 
may help to determine its origin. In this paper, we present 
observations of emission lines of CN(2-1), HCO^(3-2) and 
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C2H(3-2) towards the central galaxy and also towards a po- 
sition 20" to the east where there is strong CO emission in 
the associated filamentary structure. 

These three sp ecies were identified in a theoretical study 
IjBavet et al.ll201ll ) as tracers of regions influenced by the 
dissipation of turbulence and w aves, heating the gas and 
accelerating cosmic rays (see als o lCrawford fc F abian 1 9921 : 
iPope, Hartauist fc PittardI I2008Y Fcrlan d et all l|2009il de- 
veloped models of the optical and infrared emission filamen- 
tary regions in order to infer the rate at which mechanical 
energy is dissipated, and to determine the cosmic ra y back - 
ground produced by the interaction. iFerland et al.l (|2009l ) 
concluded that the heating rate per hydrogen nucleus due 
to dissipation or cosmic ray induced ionisation must be 300 
times higher in the optical emission filaments than in the 
local interstellar medi um. 

In their work, (.Bavet et al.ll2011h identified the species 
CN, C2H and HCO"'' as sensitive to changes in cosmic ray 
ionisation rat e or to an addit i onal source of heating such 
as dissipation. iMeiierink et al.l l)201l| ) came to a similar con- 
clusion. Here, we present the observational follow-up to that 
theoretical study; we aim to determine whether dissipation 
or cosmic ray induced ionisation dominates the heating. 

The paper is organised as follows: in Section [2] we 
present the first observations of CN(2-1), C2H(3-2) and 
HCO+(3-2) emissions in both NGC 1275 and a position 20" 
to the east where emission in the CO filamentary structure is 
known to be strong (see Figure [l]). In Section [3] we analyse 
the data and provide rough estimates of column densities 
and fractional abundances of the three molecules in NGC 
1275 and in the filamentary region. Section [4] c ompares the 
result s of the observations with the models of lBavet et al.l 
l|2011l ) to try to identify whether cosmic ray heating or dis- 
sipation is the dominant heating mechanism, or whether a 
combination of both is required. Section [S] discusses the re- 
sults and gives our conclusions. 



2 OBSERVATIONS 

The observations were performed using the James Clerk 
MaxweU Telescope (JCMT) for the detection of the CN(2-1) 
line {v = 226.874 GHz), C2H(3-2) transition (1/ = 262.004 
GHz) and HCO+(3-2) line (u = 267.557 GHz). In this first 
attempt of detection, we did not seek to resolve the fine 
structure line emission of the CN nor of the C2H but rather 
observe the emission coming from the whole (2-1) and (3-2) 
group, respectively. Every two or three hours the pointing, 
focus and calibration were performed carefully on planets 
(Mars and Jupiter) and on evolved stars. The pointing error 
was estimated to be ^ 3". 

The JCMT observations were made between August 
and December 2010 under medium weather conditions 
(t225 ~ 0.08-0.15). We employed a beam switch mode with a 
throw of 95" and used the RxA3 receiver on JCMT for de- 
tecting these three molecular emissions, coupled with the 
ACSIS digital autocorrelation spectrometer with a band- 
width of lOOOMHz because the lines were expected to be 
broad (see Figures [T] and [3]). The HPBW and the main beam 




Figure 1. Extracted from I Salome et al.l l l2008bl ). This figure 
presents the positions and size of our JCMT beams (20" grey 
circles) used for our observations relative to the integrated 
CO(2-l) intensity (white contours) seen on an Ho image from 
IConselice. Gallagher fc Wvsj l l200lh (colour scale) of NGC 1275. 
The white circle shows the Plateau de Bure Interferometer pri- 
mary beam of the Salome et al. (2008b) CO data. The black circle 
located on NGC 1275 (at (-20",-5") offsets from the centre on the 
map) corresponds to the 11.7" beam of the lR AM-30m single-dish 
telescope CO data from lSalome et al.l 1120061 ) . 

efiiciencies of the JCMT at i^ ~ 250 GHz are 20" and 0.6^, 
respectively (see the projected JCMT beam on Figure [T] as 
compared to the one corresponding to the published data 
from IRAM-30m and IRAM-Plateau de Bure Interferome- 
ter). The system temperatures ranged between 200 K and 
450 K, depending on the source (either NGC 1275 or the 
filament) and on the wavelength. The data pre-reduction 
was done using Starlink software (KAPPA, SMURF and 
STLCONVERT packages) and subsequently translated to 
CLASS format for final reduction. The reduced spectra for 
NGC 1275 and the filamentary region are seen respectively 
in Figures [2] and [3] and their resulting Gaussian fitting pa- 
rameters are displayed in Table [T] 

The C2H(3-2) line in NGC 1275 and the three lines ob- 
served in the filamentary region show a signal-to-noise ra- 
tio below the usual extragalactic cutoff of 3a. For clear de- 
tections in the filamentary region, one would require more 
on-source integration time. These observations are thus to 
be considered all as upper limits, even though in the case 
of the C2H(3-2) line in NGC 1275 and of the HCO+(3-2) 
and CN(2-1) transitions in the filament, Gaussian line pro- 
files could be fitted to the observations (see Figures [5] and 
Ol. Where fits could be obtained, the fitting parameters are 
shown in Table [T] However, note that for the lines of low 
signal-to-noise the values should be regarded, with caution, 
as indicative; these values are shown in italic in Table [T] 

To correct the integrated line intensities for beam di- 
lution effects, we have assumed optically thin emission (an- 
tennae temperature proportional to the column density in 



1 See JCMT website: |http://www.jach.hawaii.edu/JCMT/| in- 
struments/ 
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Table 1. Observational parameters and Gaussian fits parameters obtained for the data set. When the data are upper limits and we have 
not succeeded in fitting any Gaussian line profiles, we have used a symbol "— " in the table and given for the Tpeafc value a 3 it estimate 
based on the value of the observed rms. When fitting values are seen in a italic text font, it means that the observations are also to be 
considered as upper limits but it has been possible to fit a Gaussian line profile to the data (see text in Section [2]|. 





NGC 1275 


Filament 


RA(J2000) (h:m:s) 


03:19:48.20 


03:19:50.00 


DEC(J2000) (° : ' : ") 


41:30:42.0 


41:30:47.0 


Position in Figure [T] 


D 


A, B&C 


HCO+(3-2), i^ = 267.557 GHz 






Tsys (K) 


253 


315 


Integration time'^ (mins) 


20 


18 


f{Tmb dv) (K km s-i) 


18.3±0.4 


4.2±0.4 


Vpeak (km s-1) 


5178.7±5.0 


5185. 7±1 1.1 


FWHM (km s'^) 


437.5±10.0 


372.1±34.2 


Tpeafc (mK) 


39.4 


10.5 


rms (mk) 


7.9 


8.6 


CN(2-1), 1^ = 226.874 GHz 






Tsys (K) 


197 


252 


Integration time"^ (mins) 


27 


20 


/(T^b dv) (K km s-i) 


5.9±0.2 


0.7±0.2 


Vpeak (km s-1) 


5197.3±4.3 


5248. 9± 3 3. 5 


FWHM (km s'^) 


275.9±9.5 


214.2±91.0 


Tpeafc (mK) 


20.2 


3.0 


rms (mk) 


4.8 


2.8 


C2H(3-2), ly = 262.004 GHz 






Tsys (K) 


306 


442 


Integration time" (mins) 


18 


18 


f{Tmb dv) (K km s-i) 


3.1±0.2 


- 


Vpeak (km s-1) 


5126.8±7.4 


- 


FWHM (km s-i) 


180.8±n.3 


- 


Tpeafc (mK) 


16.3 


<23.4 


rms (mk) 


8.3 


7.8 



The integration time listed here is the ON-SOURCE time only. 



the upper level of the observed transition) and source sizes 
of 20" and 11" for NGC 1275 and the filament regions, re- 
spectively. These source sizes have been derived f rom the 
CO single-dis h and interferometric maps from Salome et al.l 
l|2006l . l2008al lb|) (see also Figure [T]). We have applied this 
correction before estimating the column densities and frac- 
tional abundances (but note that the data in Table [l] do not 
include any correction for beam dilution). 

We note that the HCO+(3-2) line in NGC 1275 is 
significantly broader than other lines presented here, and 
is broader than the CO(l-O) and CO(2-l) line wi d ths (o f 
about 380 k r n s~^ ) repor ted in [Bridges fc IrwinI (|l998t) : 
ISalome et"all (|2006l . 12008^ 1. The cause of this broadening 
is unclear. The range of frequencies covered by the spec- 
trum seen in Figure [2] (i.e. 267.093-268.036 GHz), includes 
only one group of lines potentially strong enough to blend 
with the HCO"'"(3-2) emission and cause this broadening: 
this is the methanol emission located around 267.5 GHz. 
However, since methanol (in any frequency range) has not 
been detected in NGC 1275, nor has it been found in any 
central galaxy of a cluster, this possible explanation remains 
hypothetical. 



Using previously published Perseus cluster molecular 
line data, we have listed in Table [5] some line ratios cor- 
rected for beam dilution effects. We have included additional 
molecular lines know n to be reliable tracer s of relatively 
high density gas from [Bridges fc IrwinI l| 19981 ) : ISalome et al.l 
t2008al ). For the filamentary regi on, we have used t he CO (3- 
2) integrated line intensity from [Bridges fc IrwinI (| 19981 ) at 
the position closest to ours ((-|-7";-|-7") as compared to 
(-|-20";+5")), without finding a better match. 



We briefiy compare our results for the Perseus Cluster 
with the results of oth er molecular line studies of galaxies. 
JGao fc Solomon! l|2004l ) give integrated line intensities (K 
km s~^) for CO in nearby galaxies ranging from a few to 
a few hundred, while for HC N the values rang e from a few 
tenths to a few tens (see also lAalto et al.ll2002f) . Our results 
(in the same units) for CN(2-1), C2H(3-2) and HCO+(3-2) 
are 5.9, 3.1 and 18.3. So, in NGC 1275 these three molec- 
ular species are typically as bright as HCN in many other 
galaxies, though not as bright as CO. 
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Table 2. Line ratios using previously publishe d data from 
[Bridges fc IrwinI l|l998h : [SaIomi et alj l|2006l . l2008ah and the ob- 
servations presented here. The line ratios have been calculated 
only after the integrated line intensities of all the observations 
used here were corrected for beam dilution effects using an as- 
sumed beam size of 20". When the symbol '-' is used, it means 
that the data are unknown. The CO(3-2) integrated line in t ensity 
used in the filamentary calculations from [Bridges fc IrwinI l|l998l ) 
does not correspond exactly to the position observed here and 
therefore has to be considered as a lower limit. 



HCO+(3-2)/CO(3-2) 

HCO+(3-2)/HCN(3-2) 

CN(2-l)/CO{3-2) 

CN(2-l)/HCN(3-2) 

C2H(3-2)/CO(3-2) 

C2H(3-2)/HCN(3-2) 



3 COLUMN DENSITY AND FRACTIONAL 
ABUNDANCE ESTIMATES 

To convert the integrated line intensities into column den- 
sities and fractional abundances, the hypothesis of a Lo- 
cal Thermodynamical Equilibrium (LTE) has been assumed 
such as, for a species X, we have: 



NGC 1275 


Filament 


2.12 


<0.53 


13.46 


- 


0.68 


<0.09 


4.31 


- 


0.36 


- 


2.28 


- 



NtotiX) = ^{X)x QiTrot) X exp 



where — '^(X) is defined by: 



Eu 



K yi 1 rot 



gu 



ix) 



L669xl0"x Z^-"'^"— 



u X (5At2) 



(1) 



(2) 



where J T,ntdVcorr corresponds to the integrated line inten- 
sities listed in Table [l] corrected from beam dilution effect, 
where Q{Trot) is the partition function at the rotational tem- 
perature Trot, where Eu is the energy of the upper level of 
the studied transition, where k is the Planck constant, and 
where S^^ is can be found, similarly to Q{Trot) values, in the 
Cologne Database for Molecular Spectroscopy (CDMqfl). 

Any more sophisticated treatment at this stage of the 
observations has been considered as meaningless since only 
one transition per molecule has been obtained per source. 
Since the current dataset of observations does not allow us 
to estimate the rotational temperature Trot i. e. rotational 
diagrams could not be constructed, we present in Table [3] 
values of the total column densities for CN, C2H and HCO''' 
as functions of Trot- The Trot values listed have been cho- 
sen to be consistent with the range of kinetic temperatures 
co mputed self- c onsist ently in the various models explored 
bv lBavet et~all (|201lf ). 

To estimate the fractional abundance of a molecule from 
its estimated total column density, Ntot{X) we have assumed 
the canonical gas:dust ratio for the Milky Way. This value 
may not be appropriate for the Perseus Cluster. We have se- 
lected Av to be equal to either 3 mag or 8 mag. This allows 
the di rect comparison with results presented in lBavet et al.l 
l|201lh . Though we cannot be certain that the material is at 
either of these visual extinctions. However, we have selected 



See the website:|http://www. astro. uni-koeln.de/cgi-bin/cdmssearch| 




-0.02 

0.03 



HCO""(3-2) 
267.557 GHz 
beam size~19.7" 



CN(2-1) 
226.874 GHz 
beam size~19.7' 




C2H(3-2) 
262.004 GHz 
beam size~19.7" 




Velocity (km/s) 

Figure 2. Molecular spectra obtained in NGC 1275, galaxy lo- 
cated at the centre of the Perseus cluster. Top panel: HC0"'"(3- 
2) line. Middle panel: CN(2-1) transition and Bottom panel: 
C2H(3-2) line. The vertical black line represents the position 
of the Visr=5125 km s"! for the CO gas in NGC 1275 (see 
Bridges & Irwin 199q). In each panel, the back curve shows the 
result of the Gaussian fit profile applied whose output parame- 
ters are listed in Table [T] (see Section [2J. Despite the fact that 
a Gaussian profile can be fitted to the C2H(3-2) spectrum, this 
observation can not be considered as a detection since its cor- 
responding signal-to-noise ratio is below the usual extragalactic 
cutoff of 3. 



them as the lower value is representative of translucent re- 
gions and the higher corresponds to a dark region that does 
not have such a high column density that its lifetime would 
necessarily be significantly limited by gravitational collapse. 
We recognize that the whole procedure that we have used is 
necessarily very crude, but it does allow rough estimates of 
fractional abundances to be made, to compare with models 
(see Section |4| . The derived column densities and fractional 
abundances, for both NGC 1275 and the filamentary region, 
are presented in Tables |3] and lU respectively. 
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HCO*(3-2) 
267.557 GHz 
beam size~19.7" 




Velocity (km/s) 

Figure 3. Molecular spectra obtained in a shifted position to 
the centre of the Perseus cluster corresponding to a filamentary 
structure. Top panel: HCO+(3-2) line, Middle panel: CN(2-1) 
transition and Bottom panel: C2H(3-2) line. See caption of Fig- 
ure (2] All the observations presented in this figures have to be 
considered as upper limits only, even if a Gaussian line profile 
can been seen for the HCO+(3-2) and the CN(2-1) lines. Indeed, 
for these three observations, the signal-to-noise ratio obtained is 
actually belowr the usual extragalactic cutoff of 3. Any Gaussian 
fit seen has thus to be considered as only indicative of a poten- 
tial detection. The velocity window seen on the spectrum of the 
C2H(3-2) line is indicative too, for marking the expected position 
of the line, if it were detected. 



COMPARISON WITH MODEL 
PREDICTIONS: CONSTRAINTS ON THE 
ORIGIN OF THE HEATING MECHANISMS 
IN THE PERSEUS CLUSTER 



We have extracted from iBavet et al.l (|201lh . the values of 
the fractional abundances of CN, C2H and HCO^ from five 
models tliat may be appropriate for reproducing the obser- 
vations presented here in both NGC 1275 and in the filament 
region. We present these values in Table O They offer five 
distinct combinations of cosmic ray ionisation rate (here- 
after (" in s^^), of additional source of heating, H, (in erg 
cm~'^ s~^) accounting for the effect of dissipation, and of 
FUV radiatio n field, /, (expressed in units of Habing - see 
iHabineHToei ). 

More precisely, Model A has high values of both H and 
( with a low value of /. Model B has low values of H , ( and 
/. Model C has low values of H and (" with a high value of 



I. Model D has a high value of H and low values of (^ and /. 
Model E has a low value oi H , a, high value of ^ and a low 
value of I (See Table [SJ . Each model explored thermal and 
chemical properties at two values of A^, i.e. Av —3 and 8 
mag. These values are intended to represent conditions near 
the edge and in the dark interior of a molecular region. 

By comparing directly the observed fractional abun- 
dances values presented in Table|4]with those modelled listed 
in Table O for the same temperature range (i.e. here, we 
compare only observed and modelled fractional abundances 
at Trot ~ Tk), we may be able to identify which heating 
process may be present and consistent with the observed 
emissions of both NGC 1275 and the filamentary region. 
A summary of our conclusions is presented in Table [6] As 
one can see, there are actually a relatively limited number 
of models which can reproduce the observations well. Here 
we arbitrarily considered that a model reproduces well an 
observation when there is a factor of 5 or less between the 
modelled and observed fractional abundances. This factor 
has been selected such as taking into account possible un- 
certainties on the detections (e.g. pointing, calibration, etc) 
as well as uncertainties in the gasidust ratio that affect the 
column densities and fractional abundances; also, the model 
results for abundances scale with the unknown metallicity. 

Inspection of Table |S] suggests that the radiation field 
has little effect on the resulting chemistry, whatever values 
are attained by H and C,. We shall ignore the FUV radiation 
field parameter in the following discussion, and focus on the 
effects of varying H and (". Note that where the observations 
supply only upper limits, we have assumed that the actual 
value is close to the upper limit. Further, we note that the 
success or failure of the models to match observations ap- 
pears to be the same for both the centre, i.e. towards NGC 
1275, and the eastern filamentary position. 

Table |6] shows that individual species may have a match 
with observation for several models. For example, CN gives 
a match to observations for Models B, E (at Av = 3 mag), 
Model D and Model C, while C2H could match the observa- 
tions for Models E and A. However, if all the molecules are 
assumed to share the same space and physical conditions, 
then we require a single model to account for all the species 
considered. Table|6]shows that there is such a model. Model 
E. It can account for CN and HCO^ at Ay — 3 mag (though 
not at 8 niag), and for C2H. The values of Trot associated 
with all three molecules is 50 K. It is encouraging that the 
same value of Trot applies to all three species. We discuss 
the implication of this finding in Section [5] 



5 DISCUSSION AND CONCLUSIONS 

The comparison in Section |4] between the molecular abun- 
dance s obtained f r om th e observations and from the mod- 
els of iBavet et al.l l|201H ) suggests that we can identify one 
model, Model E, that produces results consistent with the 
observations. This model has a low FUV radiation field, a 
low rate of heating by dissipation, but a high cosmic ray ion - 
ization rate. Unpublished data from the lBavet et al.l (|201in 
calculations show that cosmic ray heating accounts for 85% 
of all heating at either ^t, = 3 or 8 mag. 

The enhanced cosmic ray fiux may arise in_a dynami- 
cal interaction, as suggested bv lFerland et al.l l|2003 ). Those 
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authors also found that an enhanced cosmic ray ionization 
rate was required to account for the observations of the op- 
tical and infrared emitting components of the filaments. It is 
interesting, therefore, that the results from the present work 
also suggest that the regions of cluster gas probed by mil- 
limetre and sub-millimetre emissions also require a similarly 
enhanced cosmic ray ionization rate. A second inference is 
that a high heating rate from energy dissipation is not re- 
quired; in fact, as Table[6]indicates, a high heating rate from 
sources other than cosmic rays appears to inhibit a match 
between models and observations. 

In the observations, we have detected only one line in 
each of three molecular species, so a reliable analysis of the 
observational data cannot be made. It would be useful for 
the further study of galaxy cluster gas to make confirmed de- 
tections of several lines of each of the three molecular tracers 
used in this work. In the mod elling of the chem istry of the 
gas in and around NGC 1275 (|Bavet et al.ll201lh . some im- 
portant physical parameters (e.g. gas:dust ratio and metal- 
licity) are poorly known and undoubtedly have an impact on 
the model predictions. Nevertheless, the present work shows 
the value of molecular line observations of cluster gas, and it 
is interesting that the results reported here for the molecu- 
lar gas are reasonably consistent with the results for regions 
probed by the optical and infrared. The main conclusions of 
this work are that a cosmic ray ionization rate enhanced by 
a factor of at least one hundred in the cluster gas is required 
to establish the observed molecular abundances, and that 
a high heating rate from other causes appears to mitigate 
against that chemistry. In these circumstances, the 1100"*" 
emission is strong. 
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Table 3. Total column density estimates (in cm~^) obtained using the formula described in Section[3l for various values of Trot- Since 
the there are observed upper limits for both NGC 1275 and the filamentary region, the corresponding derived total column densities 
presented here are thus also to be considered as upper limits (symbol < used). The values of T^ot has been selected with respect to the 
available online values of the partition function QiTrot) and such as to be compatible with the range of kinetic temperature displayed 
by the models (see Table O. 

Trot = 2.725 K Trot = 5 K Trot = 9.375 K Trot = 18.75 K Trot = 37.5 K Trot = 50 K Trot = 75 K 

NGC 1275 

HCO+ 6.93 X 10" 1.84 x lO^^ 2.62 x 10^^ g 45 ^ 10" 6.22 x lO^^ 5.86 x lO" 5.46 x lO" 

CN 2.42 X 1014 1.91 X lO^^ 5.02 x lO^^ 2.53 x lO" 1.93 x lO^^ 1.86 x lO^^ 1.79 x lO^^ 

C2H < 9.03 X 10^5 < 1.55 X 10" < 1.74 x 10^^ < 5.40 x lO^^ < 3.23 x lO^^ < 2.94 x 10^^ < 2.66 x 10^^ 

Filament 



HCO+ < 2.10 X 10" < 5.58 x lO^^ < 7 94 ^ lo" < 2.87 x lO^^ < 1.88 x lO" < 1.78 x lO" < 1.66 x lO" 

CN < 4.68 X 10" < 3.69 x lO^^ < 9.70 x lO^^ < 4.89 x lO^^ < 3.73 x 10^^ < 3 gg x lO^^ < 3 35 ^ 10" 

C2H < 1.34 X 10^^ < 2.52 X 10" < 2.97 x 10^^ < 9.47 x lO^^ < 5.74 x lO^^ < 5.25 x 10^^ < 4.77 x 10^^ 

Table 4. Fractional abundances (with respect to the total number of hydrogen nuclei) obtained us ing hypothesis des cribed in Section [S] 
for various values of Trot and for Ay = 3 mag and 8 mag (consistently with model predictions from iBavet et alj|201ll . See text in Section 
O. See caption of Table [3] 

Trot = 2.725 K Trot = 5 K Trot = 9.375 K Trot = 18.75 K Trot = 37.5 K Trot = 50 K Trot = 75 K 

NGC 1275 



Ay = 3 


mag 
















HCO+ 




1.44 X 10"'^ 


3.84 X 10-=* 


5.45 X 10-10 


1.97 X 10-10 


1.30 X 10-10 


1.22 X 10-10 


1.14 X 10-10 


CN 




5.04 X 10-8 


3.97 X 10-s 


1.05 X 10-'-* 


5.27 X 10-10 


4.01 X 10-10 


3.87 X 10-10 


3.72 X 10-10 


C2H 




< 1.88 X 10-*^ 


< 3.23 X 10-** 


< 3.62 X 10-9 


< 1.12 X 10-9 


< 6.72 X 10-10 


< 6.12 X 10-10 


< 5.54 X 10-10 


Ay =8 


mag 
















HCO+ 




5.42 X 10"* 


1.44 X 10-^ 


2.05 X 10-10 


7.39 X 10-11 


4.86 X 10-11 


4.58 X 10-11 


4.27 X 10-11 


CN 




1.89 X 10"* 


1.49 X 10-s 


3.92 X 10-10 


1.98 X 10-10 


1.51 X 10-10 


1.45 X 10-10 


1.40 X 10-10 


C2H 




< 7.06 X 10-'' 


< 1.21 X 10"** 


< 1.36 X 10-9 


< 4.22 X 10-10 


< 2.52 X 10-10 


< 2.29 X 10-10 


< 2.08 X 10-10 


Filament 


Ay =3 


mag 
















HCO+ 




< 4.38 X 10-« 


< 1.16 X 10-3 


< 1.65 X 10-10 


< 5.97 X 10-11 


<3.93 X 10-11 


< 3.70 X 10-11 


< 3.45 X 10-11 


CN 




< 9.74 X 10-« 


< 7.68 X 10-9 


< 2.02 X 10-9 


< 1.02 X 10-9 


< 7.76 X lOio 


< 7.50 X 10-10 


< 7.00 X 10-10 


C2H 




< 2.79 X 10-6 


< 5.24 X 10-** 


<6.18 X 10-9 


< 1.97 X 10-9 


< 1.20 X 10-9 


< 1.09 X 10-9 


< 9.93 X 10-10 


Ay =8 


mag 
















HCO+ 




< 1.64 X 10-** 


< 4.36 X 10-10 


< 6.20 X 10-11 


< 2.24 X 10-11 


< 1.47 X 10-11 


< 1.39 X 10-11 


< 1.29 X 10-11 


CN 




< 3.65 X 10-« 


< 2.88 X 10-'-* 


< 7.58 X 10-10 


< 3.82 X 10-10 


< 2.91 X 10-10 


< 2.81 X 10-10 


< 2.62 X 10-10 


C2H 




< 1.05 X 10-'5 


< 1.97 X 10-** 


< 2.32 X 10-9 


< 7.40 X 10-10 


< 4.49 X 10-10 


< 4.10 X 10-10 


< 3.72 X 10-10 
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Table 5. Predicted values of the fractional abundances (with respect to the total number of hydrogen nuclei) of CN, C2H and HCO"*" 
from models seen in Bayet ct al. (2011). More precisely, Models A - E listed here correspond to Models 16, 14, 5, 18, and 12, respectively, 
in lBavet et alj 1 1201 li ) hence to models having high values of both H and Q with a low value of 7; low values of H, (^ and /; low values 
of H and C, with a high value of /; high value of H and low values of Q and / and low value of H, high value of C, and low value of /, 
respectively. See the text in Section |4] 



Model 


A„ 


Tk 




H 


C 




/ 


Frac. abun. 


Frac. abun. 


Frac. 


abun. 




(mag) 


(K) 


(crgcm- 


-3 s-1) 


(s- 


-') 




(Habing) 


CN 


C2H 


HCO+ 


A 


3 


43.8 


1.00 


X 


10-20 


5.00 X 


10- 


-15 


10 


7.33 X 10-9 


3.32 X 10-10 


4.69 X 


10-10 


A 


8 


47.5 


1.00 


X 


10-2" 


5.00 X 


10- 


-15 


10 


7.24 X 10-9 


2.74 X 10-10 


5.94 X 


10-10 


B 


3 


13.8 


1.00 


X 


10-22 


5.00 X 


10- 


-17 


10 


1.74 X 10-9 


9.01 X 10-12 


1.36 X 


10-11 


B 


8 


14.9 


1.00 


X 


10-22 


5.00 X 


10- 


-17 


10 


1.81 X 10-10 


5.22 X 10-12 


5.32 X 


10-11 


C 


3 


25.3 


1.00 


X 


10-22 


5.00 X 


10- 


-17 


100 


1.55 X 10-9 


5.50 X 10-11 


1.31 X 


10-11 


C 


8 


14.5 


1.00 


X 


10-22 


5.00 X 


10- 


-17 


100 


1.82 X 10-10 


5.28 X 10-12 


5.16 X 


10-11 


D 


3 


74.8 


1.00 


X 


10-20 


5.00 X 


10- 


-17 


10 


2.51 X 10-9 


8.82 X 10-11 


1.65 X 


10-10 


D 


8 


63.3 


1.00 


X 


10-20 


5.00 X 


10- 


-17 


10 


8.61 X 10-10 


5.43 X 10-12 


2.62 X 


10-10 


E 


3 


51.0 


1.00 


X 


10-22 


5.00 X 


10- 


-15 


10 


1.82 X 10-9 


2.68 X 10-10 


1.37 X 


10-10 


E 


8 


60.2 


1.00 


X 


10-22 


5.00 X 


10- 


-15 


10 


1.99 X 10-9 


2.95 X 10-10 


2.94 X 


10-10 



Table 6. Summary of the conclusions obtained when comparing the estimates of the fractional abundances derived from the observations 
to the model predictions (see Section |4j. We have assumed that models and observations are in agreement (symbol '+' used) when their 
values differ from less than or equal to a factor of 5 difference. Otherwise, the symbol '-' is used. '(3)' or '(8)' means that the fractional 
abundances at Ay = 3 mag or at Ay = 8 mag are well reproduced by the models (but not at both Ays), respectively. 



Model 


C2H 


C2H 


CN 


CN 


HCO+ 


HCO+ 


Parameters 


centre 


filament 


centre 


filament 


centre 


filament 


Model A: I low, H high, ( high 


+ 


+ 


- 


- 


(3) 


- 


Model B: / low, H low, C, low 


- 


- 


+ 


+ 


(8) 


(8) 


Model C: / high, H low, ( low 


- 


- 


+ 


+ 


(8) 


(8) 


Model D: / low, H high, ( low 


- 


- 


+ 


+ 


(3) 


(3) 


Model E: / low, H low, f high 


+ 


+ 


(3) 


(3) 


(3) 


(3) 



